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Georgi  Zlatanov 


Introduction 


In  recent  years  a rapid  development  of  science  and  tech- 
nology is  observable.  In  the  field  of  geodesy  this  develop- 
ment is  characterized  by: 

- The  introduction  and  rapid  spread  of  electromagnetic 
stadia 

- The  individuation  of  a new  geodetic  field  - cosmic 
geodesy 

- The  versatile  implementation  of  electronic  technology 
to  geodetic  techniques 

These  new  elements  in  geodesic  theory  and  practice  have 
enabled  new  methods  for  the  reciprocal  establishing  of  the 
location  of  points  on  the  surface  of  the  Earth.  These  methods 
are  essentially  different  from  classical  methods. 

In  classical  geodesy  all  dimensions  were  divided  into 
two  basic  groups:  horizontal  dimensions  and  vertical  dimen- 
sions. Classical  geodesy  required  a reference  surface,  upon 
which  were  reduced  all  measurement  results. 

The  problem  of  measurement  reduction  was  one  of  the  basic 
problems  of  classical  geodesy.  The  accuracy  with  which  reduc- 
tion corrections  could  be  determined  constituted,  in  many 
cases,  a barrier  determining  the  accuracy  of  the  geodesic  net. 

Modern  technology  allowed  that  barrier  to  be  overcome. 

An  important  role  in  this  was  played  by  satellite  geodesy, 
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which  in  a short  time  achieved  a high  level  of  development. 

An  essential  share  in  this  development  was  contributed  by 
Polish  geodesists. 

As  is  well  known,  cosmic  triangulation  can  be  thick- 
ened, with  the  aid  of  the  so-called  balloon  triangulation. 
As  a result  of  this  thickening  the  geocentric  coordinates  of 
a large  number  of  points  on  the  Earth's  surface  can  be 
gained.  Consequently,  it  is  entirely  possible  to  obtain  the 
geocentric  coordinates  of  points  lying  20C-300  km  apart, 
v/ithout  the  necessity  of  reduction  (in  the  classical  geodesic 
sense)  results  of  direct  measurement,  and  this  without  using 
a reference  surface.  These  new  possibilities,  essentially 
differing  from  classical  means,  require  the  discovery  of  new 
methods  for  establishing  basic  geodesic  nets. 

In  this  current  wo rk  a method  for  the  mathematical  work- 
ing out  of  the  geodetic  net  is  considered,  based  on  three- 
dimensional  geodesy.  Adjustment  and  calculation  of  the  net 
is  performed  on  a three-dimensional  Cartesion  coordinate 
system,  and  not  on  a reference  surface.  By  this  method  the 
difficulties  resulting  from  the  reduction  problem  are  avoided. 
There  arises,  however,  the  necessity  of  utilizing  physical 
reduction,  and  the  especial  need  to  consider  the  effect  of 
refraction. 

Due  to  the  character  of  the  work,  it  is  assumed,  that 
the  angles  of  refraction  are  determined  beforehand  (e.g., 
with  the  aid  of  a refractometer) , and  thus,  that  the  zenith 
distances  are  free  from  the  influence  of  refraction. 


In  three-dimensional  geodesy  it's  necessary  to  make  an 
aggregate  adjustment  for  the  effects  of  the  measurement  of 
the  vertical  and  horizontal  angles.  There  thus  arises  the 
possibility  of  making  calculations  in  a three-dimensional, 
Cartesian  coordinate  system.  Ko  difficulties  are  presented, 
in  this  case,  including  the  effects  of  length  measurements. 

The  idea  of  three-dimensional  non-reduced  geodesy  arose 
in  the  middle  of  the  last  century  in  the  work  of  Willarseau 
and  Bruns.  At  that  time,  however,  its  development  could  not 
be  continued.  Only  irf  the  1950 's  do  we  encounter  the  first 
systematic  investigations  dedicated  to  three-dimensional 
geodetic.  The  groundwork  for  a geometric  approach  to  the 
question  of  determining  the  figure  of  the  Earth  was  laid  by 
K.  S.  Fo^odienski.  A series  of  basic  mathematical  equations 
for  three-dimensional  geodesy  can  be  found  in  the  work  of 
Botin  and  Dufor. 

After  formulating  the  basic  laws  of  cosmic  geodesy,  and 
after  working  out  the  geometric  methods  for  investigating  the 
Earth's  figure,  there  began  researches  regarding  problems 
arising  in  conducting  local  spatial  triangulation. 

The  set  of  problems  included  here  can  be  divided  into 
two  groups: 

1.  The  effects  of  refraction  on  the  results  of  zenith 
distance  measurements 

2.  Methods  of  adjustment  for  spatial  networks 

Whereas  the  question  of  the  effect  of  refraction  finds 
a diverse  and  occasionally  contradictory  variety  of  inter- 
pretations, methods  for  adjustment  of  spatial  networks  can 
be  assigned  to  two  groups: 
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1.  Separate  adjustment  of  the  vertical  and  horizontal 
measurement  results 

2.  Aggregate,  i.e.,  three-dimensional  adjustment  of 
vertical  and  horizontal  measurement  results. 

By  separately  adjusting  these  first  is  determined  the 
horizontal  coordinates  of  the  network  points,  and  next  the 
third  coordinate,  together  with  the  unknown  perpendicular 
deviation.  The  means  of  adjustment,  belonging  to  this  group, 
do  not  differ  essentially  from  the  classical  geodetic  methods 
discussed  in  textbooks.  They  represent  a single  minor  modi- 
fication of  the  known  methods.  The  essential  difference  in 
relation  to  classical  means,  is  shown  by  the  methods  for  estab- 
lishing aggregate  adjustment  of  the  horizontal  and  vertical 
angle  measurement  results.  In  this  adjustment  perpendicular 
deviations  also  play  a role.  It  results  that  measurements 
carried  out  on  the  surface  of  the  Earth  are  directly  connected 
with  the  direction  of  the  perpendicular  line  at  the  observa- 
tion stations.  There  arises  thus  the  necessity  to  reduce 
them  to  a uniform  system  of  coordinates. 

In  the  current  work  two  algorithms  for  spatial  triangu- 
lation adjustment  are  presented;  in  the  first  of  these  multi- 
group adjustments  by  the  conditional  method  is  used  with 
unknowns,  and  in  the  second  — the  parameter  methoc.  of  adjust- 
ment for  correlated  quantities. 
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Whereas  in  cosmic  geodesy  the  axes  of  coordinates  of  a 

topocentric  system  of  coordinates  is  parallel  to  the 

axis  of  a geocentric  system,  in  the  case  under  examination  — 

topocentric  systems  of  coordinates  are  used.  Closely  connected 

with  the  positions  of  the  theodolite  at  the  time  of  observation. 

, y , z y 

Thus,  for  example,  the  topocentric  coordinate  system  (X,  Y,  -2)^. 
(drawing  1)  can  be  defined  as  follows:  axis  Z is  directed 
towards  the  zenith  (along  the  horizontal  axis  of  the  instru- 
ment) ; oxis  X lies  in  a horizontal  plane  and  is  directed 
towards  a zero  indicator  of  the  horizontal  circle;  and  axis 
% is  turned  90°  in  relation  to  axis  i*,  in  a clockwise  direc- 
tion. Such  a topocentric  coordinate  system  is  a left-handed 
system. 


-/ 


R..,,  - horizontal  angle 

between  axis  .iH'and 
direction  M-N 

P.™  - zenith  distance  of 
""  direction  K-N 

Dkii  - distance  F-N 


The  position  of  the  arbitrary  point  N in  the  topocen- 
tric coordinate  system  for  station  K can  be  defined  by 
polar  coordinates  (R,  , D)  and  by  the  rectangular  coordi- 
nates (-%,*  ’ £-)w . 

The  relationship  between  these  coordinates  is  as  follows: 

■ DioCoi/Ik,  ^ | | ^ 

sin/?**; 

zV'  m nm„  C 
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Prom  (I  1 ) we  can  obtain  the  reciprocal  relationship: 


tg^M*  = 


For  the  direction  cosines  MN  in  a local  coordinate  system 
we  have: 

".mv  — ■ = cos  AJjmv  * sin 

^w.V 


}.<jn 

f>MN  = --  — = sin ' sin^.MvI 

D\in 

,(*/> 

f»«  = — — = COS0M*. 


C'  •*> 


If  we  designate  the  versors  of  the  topocentric  system 
of  coordinates  by  (e^,  el,,  then  for  the  unit  direc- 

tion vector  MN,  after  using  the  general  cracovian  calculus 
operators,  we  can  write: 

»««  = aMV<’’i*,,  + hMV<JW>  + CMvijMl  = «MW  ' t(«,  (.  I • 


Vector  r^  can  also  be  presented  in  a geocentric  coor- 
dinate system  (X,  Y,  Z).  If  the  versors  of  this  system  are 
designated  by  E^,  E2,  E^,  then  rj^  can  be  written  as  follows: 


where: 


run  = Am*  • E,+Bum  ’ ^i  + Cm*  • £,  = A M*  • tE, 


X*-X« 


Om*  * 


_ Zm-2m 

C«*  = — 

Dun 


0*1 


O.tl 


/)«»=  + /(**  - a'm)1 + ( r*  - yiF+  (Z*  - Zm)*; 
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From  (1.4-)  and  (1.5)  it  results: 


f ns  — tins  * re  A,  =»  Ajka  • tE. 


OD 


We  accept  that  the  relation  between  e^fl/wE  is  repre- 
sented in  the  following  form: 

0*1 


E = u*  • Ttu. 
C3-3) 


where  is  an  orthogonal  cracovian.  It  is  designated  by 
three  topocentric  orientation  angles  of  the  topocentric 
coordinate  system,  namely: 


qp(M) 

A(m) 


- the  astronomical  azimat  of  axis  (of  a direc- 
tion corresponding  to  the  zero  indicator  on  the 
horizontal  circle  of  the  theodolite  at  point  M) , 

- astronomical  width  at  point  FI , 

- astronomical  length  at  point  M , 


For  elements  Uj,  the  following  relations  result: 

"n  “ -cosa0s*iiycos J-sina0sinil, 

Hi ■ =*  -cosa0sinysin  A + sinaoCOS^., 
iii,  --  cos  au  cosy, 

«u  =•  sina0sinycos  i — cOs«nsin  A , 

««  c sin«„sinysin A + cosaoCos A, 

«ji  - -sina0cosy; 

«u  = cosycos^, 

«ji  * cosy  sin  X , 

«jj  = siny. 


In  these  formulae  the  factor  M is  disregarded. 
From  (1.7)  and  (1.8)  it  results: 


f 
1 1 

I 


*M.V  — * tUjW  I 

(13)  (1.3)  (3.3) 

A MN  = dh4N  * UM 
(13)  (1.3)  (3.3) 


HA&) 

0-‘Oa) 


The  fr-~m  of  the  developed  formula  (1.10)  appears  as  follows: 


- 'I»n'<\'?)+B>ikII,">  + CunU<“\ 

*MI»  » AMN‘<\l>+BMNUii,f)+CttltU}i\ 


CI»^ 


The  relationship  (1.11)  presents  these  conditional 
equations  between  the  results  of  observations  R^j- , /3f.T,  , and 
coordinates  (X,  Y,  Z)MM , as  well  as  the  orientation  angles 
C°<0,  Between  these  three  equations  there  a. theses  the 

' ^ (( frwt ^ squaW.  on)-}  a**+l>l,M+c\tH  -.1,  ail  wo  oan  only 

two  of  theJ^v'f^or  example,  the  terms  and  ca»n  6c  inosrp^atcJ  m 

Before  approaching  adjustment,  v/e  must  present  rela- 
tionship (1.11)  in  linear  form. 

let  us  use  these  notations: 


Bum  — Rmn  + • Pun  * Pmn+  h<A(\ 

/ / / 

OMN  — COS  RMN&inpMH, 

hMN  = sin/tM/vsin/?MW> 

= co$Pun; 

a»o  = *«o  + . y'.w  » Vm  + d<fu , AM  = Am  + d).„ 

x*  =■  xl+dx*.  yK  = yZ+<jyk,  zK  m z£-mz* 


(A.'=A/,Arj 


Amn 


xn-x^m 

n°N,N 


flMW 


Dun 


Cun 


Zi-Zl, 

0°UN 


O(M)  _ , (M>,  0 0 ,0  . 

uij  = utj  (a,f0, 9Pv,  Am). 


(mu) 

(l.ixk) 


In  the  above  relationships  the  indicator  "prim"  desig- 
nates measured  values,  and  indicator  "0"  — approximate  values. 
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The  linear  form  of  equations  (1.11)  is  as  follov/s: 


''"A/.v  / „ <‘Xh,s  „ ('Bun  „ cCUN\ 

TrVs  ' MN  ' +"2‘  iX„  ,U  BXu)  ( ‘,  ,,)  + 

- (">> y;v  f "5 1 + «Si  <WW>W  - Lf, + „},  fgs)  («„ - «,) + 

' t,'  .w  okm  / \ tz„  oz„  ax u I 

-(sit,*:  1 ;°  ^"il  I /■'<’  d"*1  \ I («>  / ,0  ^"u  ..0  ^“11  ,.o  ^u><  \ . <A»>  , 

' MN^’‘  mn^+Ca,%«<s>)  ^ v*'*)  + 

~ (■'*"'  + a-“'w  ^(Af, + c*"  £tS))  d/<M)  f “ °-  (,-,3s 


£/>MAT  ,,  / „ t/f*. 

■ ' *'•*»  + g ' h0w.,-(«?J-T  + *J,  T.„  - - 

(Kun  ' Pun  \ (Xu  Lx„  cX„ 

- * *4.  -Si  * <4. (•<&.£ 

(.o  /wii  ro  o rMji  \ 

v*  <•;.<«> + *w\-  a..w, + c-*" k(*'.h'*°+  *'"*  “ 0 


* . n BBUs  . ('Cun 

eX  + ,<JJ  'i  V +"”j  \ - 

<A»«  i** 


6"»i  ,n  ™J3 
< /.'*<, + 'M,v  Afj 


' 1 ^ 4 <».£ 

«<'">  + "HT”'.  + C',v k(«.)'/;,M,+  M“*  “ 


iCM»  \ 

^~}Wu-BXk)  + 

« / 

*21  ,.0  £?«j  j \ 

(1.13b) 


" 5«>  * fl”v  +c«. ('ll. 

' (A*  Sm.  + at»  *i“w,  + C“'"  rfA<Ml+^ 


‘fMV  ./  / o „ f'n«v  „ firw*\ 

:g  ("f»  J~y  +"23  +»31 

' \ 'Am  Mm  eJf„  / 

C'w  O’u  / \ tvf«  t/Z,, 


, — -lwz„-<ttv()  + 

' ^ ',w\/  rv  \ I a ( ^mn  . (‘Bun  . &Cun\ 

Trm  )<■"»-”')-(«*  fz-M  +‘‘°»  «5jr+"*»  + 

--  +cZ,„ 

i«,  / \ 1^“)  C-y'M)  T 

n rwij  \ , 

f'M" ^fio)  ^"v  - °-  (1.13c) 


AIA  =•  0. 


I 


The  partial  derivatives  occurinp-  in  equation  (1.1J) 
are  expressed  as  follows: 


t'U.MN  . , 

' , — — —SIM 

< Khx 


C’bsiH  • 

, «=  cos/fMjvsin/J«v. 

cjKmh 


('ft  SIN 

-sin  ft's,  N. 


cos  Nun  cos  Pun  , 

siiiKM.vcos/JM*; 


< -*WV 

_ .WS#*4  I/.WV 

1 

1 Yain  ‘ sIXmn 

dC«N 

<A'« 

<v 

' v„ 

Pun 

Ww  =' 

Ohn 

< A „V 

»^°/.v  U?,.v 

/)(>J 

t'MN 

• W>i.v 

< > w 

\Xmh+  1?S is 
PlK1S 

^c.,„ 

S3 

i Yu 

AZ%„  AYin 
»Vn 

< A m M 

< «,,! 

JYiur  \/us  . 

<'CMN 

Axttn+  1 >'w  v 

<A, 

/)<»» 

D US 

<A, 

Pus 

^ u 

»M/V 

<"i  x 

< «0 

l,t“  « _„o 

..  "I2» 

fo« 

Pits  1 

= 0; 
5«0 

'"|2 

( r{, 

- t«?J. 

CHll 

- ' * ~"ll. 
cao 

- 0; 
c'a0 

< U |3 

* a(> 

II 

+ 

*So 

itti  i 

- «=  -«,3, 

COLq 

n 

-- — = 0 

£oo 

-II,  ,COs/.n  , 


-«ijCosA°, 


= -i/2jsin/.". 


— = -ttjjCOS/." 

(sty 


~ Mjjsin/.0; 


— //jj  cos  a0, 


-«r  - 


,(,«K  “ nM*~  A'sihii'},- B°in»21 
H'mn  = b\ift- /ihsii^t-  Bmniiii 
W'\Vv  = CM* /v«?3 


-Chsiii,, 

-citNii'u. 
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Utilizing  the  aid  of  EMC , the  numerical  values  of  these 
derivatives  can  be  easily  determined  by  the  i.umerical  differ- 
entiation method.  In  the  case  of  the  zenith  distance  not 
being  freed  from  the  effect  of  refraction  beforehand,  it  is 
necessary  to  include  in  the  adjustment  the  coefficient  of 

JL 

refraction, JKyj  as  an  unknown.  To  this  end  in  equation  (1.13) 
it  is  necessary  to  use  a substitution: 


~ h>Mjv  + 


1 Hu* 

2 V?ir 


(1.15) 


correction  for  measured  value  r\*l 
refraction  coefficient  for  position  K 
mean  radius  of  the  Earth 

As  we  already  know,  for  each  direction  MN  it  is  necessary 
to  employ  equations,  using  two  among  the  three  relationships 
(1.13).  Using  the  conditional  method  for  multi-group  adjust- 
ment with  unknowns  (8) , this  pair  of  equations  will  be  treated 
as  a separate  group. 


Where : 


V 


Pm 

K-ai 

far 


From  each  pair  of  conditional  equations  two  are  desig- 
nated for  equivalent  corrections,  relating  to  [vJok.d  vjnkmowkj-] 
and  independent  fictitious  measurements.  These  equations 
then  undergo  a general  operation,  as  a result  of  which  all 
unknowns  (X,  Y,  Z «x0)^  are  obtained,  also  the  cor- 

rections j Vp  for  each  n easured  direction. 

We  consider  that  including  in  the  adjustment  <f>  X 
for  each  position,  increases  to  a large  degree  the  dimensions 
of  the  system  of  normal  equations, and  therefore  it  would  be 
worthwhile  to  eliminate  these  unknowns  beforehand. 
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II 


A problem  for  geodetic  nets  installed  for  use  in  hydrau- 
lic construction  in  mountainous  regions,  or  in  the  construc- 
tion of  long  tunnels  is  assuring  high  accuracy  in  the  corre- 
lated position  of  a network  of  points. 

i 

Using  the  classical  methods  of  determining  the  effects 
of  the  measurement  of  angles  and  lengths  it  is  essential  to 
know  exactly  the  deviations  from  the  perpendicular  and  the 

geoid  interstices  from  the  elipsoid. 

I •. 

i 

Introduction  to  geodetic  calculations  of  modern  concepts 
allows  the  analytic  presentations  of  the  nets  elements  in 
three-dimensional  Cartesian  space.  Thanks  to  this,  the  neces- 
sity of  including  corrections  for  measured  quantities  in  order 

I to  transfer  to  a reference  surface,  is  eliminated. 

For  the  purpose  of  eliminating  the  effect  of  errors  in 

the  initial  data,  nets  of  high  precision,  the  rules  are 

> 

regarded  as  independent,  i . e . , not  connected  with  given  points. 

j 

The  adjustment  algorithm  considered  in  Fart  I requires 
the  existence  of  at  least  three  connected  points,  and  thus 
can  not  be  used  for  adjustment  of  an  independent  net. 

i 

Another  algorithm  is  presented  below,  which  can  be  used 
equally  in  the  consolidation  of  spatial  nets,  and  in  the 
adjustment  of  independent  nets. 

■ 

In  order  to  adjust  local  spatial  nets,  a basic  problem  is 
designating  suitable  initial  data,  which  would  determine  the 
orientation  elements  of  the  net  without  its  distortion.  There  . 
are  six  orientation  elements.  They  are  the  coordinates  of  the 
' . basic  point  (ffi  ^4—*-) 0 and  the  three  Euler  angles.  To  these 

six  orientation  elements  we  must  add  the  net  scale. 

I -!2- 

1 

\ 


The  above  seven  parameters  can  be  replaced  b;  seven  other 
independent  quantities.  There  can  be,  for  example,  the  rect- 
angular coordinates  of  two  points  and  one  coordinate  of  the 
third  point.  With  the  existence  of  three  data  points,  two 
supernumery  elements  then  appear,  which  cause  a net  deforma- 
tion. 

The  designation  of  only  one  geocentric  coordinate  for  the 
net  point  in  order  to  eliminate  its  deformation,  presents 
great  difficulties.  The  conclusion  results,  that  the  means 
of  adjustment,  well  known  in  the  literature,  which  uses  geo- 
centric coordinates  (X,  Y,  Z)  cannot  be  used  in  relation  to 
unconnected  spatial  nets. 

For  adjusting  local  spatial  nets  the  conditional  method 
can  be  used,  but  great  difficulty  arises  in  this  regard  in 
setting-up  conditional  equations.  In  regard  to  this,  in 
this  work  we  have  concentrated  on  the  parameter  method.  The 
unknowns  (X,  Y,  Z)  are  converted  in  advance  to  the  geodetic 
coordinates  (B,  L.  H) . 

Adjusting  the  net  begins  by  setting  up  the  equations 
expressing  certain  quantities,  independent  of  the  orienta- 
tion station.  The  first  time  — as  a function  of  geocentric 
coordinate  stations  and  observed  net  point,  and  the  second 
time  — as  a function  of  horizontal  and  vertical  angles, 
measured  at  the  same  station.  The  functions  spoken  of 
should  be  linearly  independent  and  should  simultaneously 
define  the  position  of  the  directions  measured  at  a given 
station. 

In  order  to  make  the  selection  of  suitable  functions, 
let  us  consider  the  directions  from  station  M to  observed 
points  I,  J,  and  K. 
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Let  us  imagine,  that  point  M is  the  center  of  a sphere 
of  elementary  radius.  Points  i,  j,  and  k (Drawing  2)  are 
then  points  intersecting  directions  MI,  MJ,  and  FK  from  the 
surface  of  the  sphere.  The  straight  sections  in  Drawing  2 
are  in  reality  arcs  of  large  circles,  corresponding  to  spa- 
tial angles.  Bearing  in  mind  that  points  M,  I,  J,  and  K are 
points  found  on  the  Earth's  surface,  the  spherical  triangle 
(i,  j,  k)  is  very  oblate  (the  degree  of  angle  at  point  K is 
nearly  180°).  In  this  situation  the  spherical  distances 
(i,j),  (i,k),  and  (k,j)  are  almost  linearly  dependent,  or 

, (2.1) 


and  thus  they  can  not  reliably  determine  the  position  of 
directions  MI , MJ  , and  MK. 


Drawing  2 


It  also  results  that  from  three  noted  spherical  distances 

for  adjustment,  only  two  are  used,  for  example,  (i,k)  and 

(k,j).  These  quantities  complement  the  spherical  altitude 
ki 

leading  from  point  R. 
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It  also  results  that  from  the  three  noted  spherical  dis- 
tances for  adjustment,  only  two  are  used,  for  example,  (i,k) 

and  (k,j).  These  quantities  complement  the  spherical  altitude 
ki 

L^.  leading  from  point  R . 

¥ 

Let  us  designate 


(tit)  * COJ  (/.  k), 

Utj)  - co »(*,/), 

(h*j)  =■  sin(/fy). 


(2.2) 


These  quantities  (2.2)  can  be  expressed  by  the  elementary 
vectors  r^,  "r^,  and  rk,  for  directions  MI,  M J , and  MK. 


Then  we  have 


(#«) 

(St,) 

(tij) 


Cd,  d). 
(r't.f}). 
(rk  r,.  rj) 
l(d  x ft)  | 


(2.3) 


For  quantities  (2.3)  the  following  dependencies  are  appropriate 


(■*«)  “ (*tl)i 

(Stj)  = (s,t),  (2.4) 

(ktj)  - -<*,')■ 


Using  quantities  of  type  (s)  and  (h)  as  adjustment  quan- 
tities does  not  create  any  essential  difficulties.  The  im- 
portant concern  here  is  the  correct  designation  of  the  corre- 
lation cracovian  of  these  quantities,  and  considering  this 
cracovian  during  the  adjustment  itself.  (6) 
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The  number  of  quantities  of  type  (s)  and  (h)  for  examin- 
ing station  M,  in  an  example  when  the  horizontal  directions 
are  measured,  and  the  zenith  angles  to  points  ^X^v^lII  equal 
(2k  - *),  in  which  (k  - 1)  ~ of  type  (s)  and  (k  - 2)  — of 
type  (h). 

In  setting-up  the  correction  equations  for  quantities 
of  type  (s)  and  (h)  formulae  are  used  which  present  elementary 
vectors  r^,  r ^ , and  r^  in  a local,  (1.4)  and  geocentric  (1.5) 
system. 

Keeping  the  local  system  in  mind,  we  have: 


(/<?>>,. 


fh  (pi.  \ Rh  R i ' W> 


(2.5a) 

(2.5b) 


and  for  the  geocentric  system: 


YM,Z%>\X,.Y„Z,\Xk,Yk,  Zt)  - 

- Hu\B„L„  H,;B, . 
Y„,Z»;X,.  Y,.Z„X„  Yj.Zj.Xl,  Y..ZJ  - 

- / «.  l.„  //,;  B„  l.j,  Hj,  B„,  U.  //«). 


(2.6a) 

(2.6b) 


The  geocentric  coordinates  (X,  Y,  Z)p  and  geodetic 
(B,  L,  H)p  for  P - I,  J,  K,  M,  are  connected  with  the  well- 
known  formula : 


X,  - (V,  + cos cost., 

Yr  - H,)cottt,tinLr, 

X,  = (/V,  t H,) tin tir- tlN,tinBr, 


(2.7) 


where  N 


P 


e(l  - e2ein2Bp)-1/2 


is  at  the  radius  of  curvature 


of  a transverse  section,  e — is  the  eccentric  of  the  mere- 
dionnl  elipse. 
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Functions  (?.r-)  and  (2.6)  are  obtained  as  a result  of 
f2.J)  and  by  considering  (1.5N  and  (1.6N.  The  final  form 
of  the  functional  dependencies  and  in  which  appear 

thegeodetic  coordinates  (B,  L,  H) , points  (M,  I,  J,  K)  is 
determined  resulting  from  (2.3),  (1.5) » and  (1.6),  and  con- 
verting in  (1.6)  coordinates  (X,  Y,  Z)  to  geodetic  (B,  L,  H) 
with  the  aid  of  relations  (2.7). 

In  order  to  obtain  the  correction  equation,  these  sub- 
stitutions are  used: 


H K-v+Vii 

X = X°  + dX,  Y-Y°+dY.  Z-  Z'+dZ; 

B » H°  + dB,  /-  - L°  + dL,  H » H°+dH ; 


(2.8a) 

(2.8b) 

(2.8c) 


and  (2.5a)  and  (2.6)  develop  in  the  Taylor  series,  limited  to 
linear  expressions. 


Then  we  have 


(*m)i«  =*  (*«)'+  Y(j„), 
I'll - OitjY+VOiif). 


(2.9a) 

(2.9b) 


where 


(*«)'  = rUti.p'k  .ft'i.Kk). 

</'*>)'»  ft),  Rj.JtO. 


(2.10a) 

(2.10b) 


V(Sik)  and  V(h^j)  are  prospective  corrections  for  quantities 
(Sik)  and  (J^)  (functions  of  directly  measured  quantities, 
subject  to  adjustment),  whose  differential  forms  are  as  fol- 
lows: 


. t (■«.»)  H*. «)  „ t 

^ “ i*rK*-+  ™r^*+ -*r  K'<+ is~  '** 

,(/.?>)  „ , wfy)  „ ^W'ifLv  ..  . wto 


(2.11a) 


, Wijt  v , Wij)  c'n,'Lvm  4._  ‘i'-l'*  K*.  f o liv'i 
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..  ■ — . 


From  the  other  side: 


(2.12a) 

(2.12b) 


I 


| 


Where: 

C2.I3<0 

C2.»3W> 

\ 


i 

1. 

Partial  derivatives  appearing  in  formulae  (2.11)  and 
(2. 14)  are  obtained  on  the  basis  of  formulae  (2.6a)  and  (2.6b). 

Using  the  dependencies 

fa**))*  » (in)..., 

or  (2.9)  and  (2.12)  we  obtain 

y(sa)  -Mu, .)+(*&>•-(*«)'. 

►'(Ay)  - KAy)+(Ay)“-(Ay)'. 
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(2.15a) 

(2.15b) 


(2.16a) 

(2.16b) 


( V,.)0  - Yfr,  zh  ■ X?,  Yl,  Z?;  x£,  Y[.  Z[)  = 

- i-°«,  //Sr.  B?,  /.?,  //P;  Bl,  Lt  Z£), 

(*y)°  = A'Ut.  Yl^Z'h  -.Xi,  Yf.Zf.Xj,  Yf.Zf  -Xt,  Yl  Zt)  - 
= '/>?<//&,  /-Sr.  //Sr  i «!’.  /.P.  //“;  Bf,  Lj,  Hi,  Bl.  L°k,  H°k). 

Us, .)  <(in)  «<*n)  /'(in)  /(in)  /(in) 

"XM  + ...  ")%t+  . d7.k,+  dX,+  i!Yi+— dZ,+ 

iXs,  tY„  (/.„  ex,  BY,  bZ, 

■ IV  x.  < II-  . ^*u)  /'(ill)  /(i,l>  /(in) 

oXk  fY,  cZk  iB „ /t*  c//„ 

/(in)  . <(*.»)  /(in)  / (in)  /(in)  ,,  t(i.») 

+v»,  ■*+ 

'4- ‘!l  gu.qiu. 

f'M  p/m  C'A4  £}| 


t (Ay) 


< (Ay) 
/>/ 


-i  ~r..' ' «/*o + ~^y-  r/zy + ^ 


*^rfjr.+  ‘$LdYk+»“bLdZk 

Cslji  C lj  i 


♦ ' ,'ir  * *£*  ♦ /&». 


i 


Expressions  (2.16a)  and  (2.16b)  are  really  equations 
for  correcting  quantities  (s)  and  (h) . In  order  to  solve 
these  equations,  as  already  noted,  it  is  necessary  to  use 
t;  e principle  of  adjustment  for  dependent  quantities  by  the 
indirect  (parametric)  method  (6). 


There  is,  however,  another  possible  approach.  Using 
formulae  (2.11a)  and  (2.11b)  the  quantities  V(s.,  ) and 

y XJC 

V(h/j)  appearing  in  equations  (2.16a)  and  (2.16b)  can  be 

replaced  by  corrections  and  ^(A)  to  directly  measured 

quantities . 

i 

Adjusting  a spatial  net  in  this  case  leads  to  another 
standard  problem  of  the  method  of  smallest  squares  — of 
conditional  adjustment  with  unknowns. 


Depending  on  the  means  of  presentation  (S.v)  and  (h^.) 

IK  1 J 

in  adjusting,  appear  as  cartesian  unknown  geocentric  coor- 
dinates (X,  Y,  Z),  or  geodetic  coordinates  (B,  L,  H)  of  net 
points. 


Installing  a spatial  net  of  high  accuracy  — in  order 
to  eliminate  distortion  of  the  net  produced  by  errors  in 
initial  data  — we  accept  as  initial  data  the  geodetic  coor- 
dinates (B,  L,  H)  of  two  points  and  the  altitude^of  the 
third  point. 


Independent  of  what  was  presented  above,  a practical 
solution  of  the  problem  encounters  difficulties  of  a prac- 
tical kind,  connected  with  the  designation  of  the  coordinates 
(B,  L,  H)  of  two  initial  points,  and  doing  so  in  such  a 
manner  that  they  would  correspond  accurately  to  the  designated 
length  between  these  two  points.  This  is  a problem  which 
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deserves  especially  close  attention.  Here  we  shall  limit 
ourselves  only  to  a discussion  of  one  of  the  substantive 
cases  — from  a practical  point  of  view. 

Let  there  be,  for  one  of  the  net  points,  which  we 
shall  call  the  beginning,  geodetic  coordinates  (B,  L,  I1)Q. 
Let  us  assume  that  we  possess  the  geodetic  azimuth  A q-^,  the 
zenith  distance  ^Q-p  and  the  incline  distance  to 
adjacent  point  1.  The  question  arises,  how  can  the  simple 
geodetic  problem  be  solved  in  space  — designating  the  geo- 
detic coordinates  (X,  Y,  Z)j? 

Here  we  shall  utilize  the  formula: 


'S^o 


cos U(Y,  - >'„)-sin£o(*i  -X») 


cos  fl0(Z,  — Z.)  ~ sin  So  [cos  A(^i  — A*o)+  sinlot^i 


^ ^ cos  R0[coi  Lo{X , — A"o)  + sin  Lo(Yj  — Yo)]+iinB0(Zl  — Zo) 

cospoi  - ]/{Xl  - Xo)1  + Zo)1  : 

Dl  ,=(*.-  Xo)1  + (Y,—  Yo)1  + (Z,  - Zo)1 


(2.17) 


In  the  above  formulae  we  can  regard  as  known  quantities 
(X,  Y,  Z)Q,  (B,  L,  H)q,  A01,  D01,  and  as  unknowns  — 

(X,  Y,  Z). 

In  other  words,  we  have  three  nonlinear  equations  with 
throe  unknowns.  Solving  them,  for  example  by  the  Newton 
method,  there  is  obtained  the  required  exactitude  (X,  Y,  Z)^ 
Next  with  the  aid  of  formula  (2.7)  the  transformation 
(X,  Y,  Z)(  (B,  L,  H)j  is  performed,  thanks  to  which  the 

problem  having  as  a goal  the  determination  of  coordinates 
(B,  L,  H)  for  two  points,  in  order  to  determine  sufficiently 
accurately  the  net  scale,  is  solved  from  both  the  practical 
and  theoretical  points  of  view. 
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By  including  also  the  ellipsoidal  height  (II)  of  the 
t ird  point,  we  obtain  the  parameters  essential  for  the 
entire  net. 


In  the  case  of  our  having  other  astronomic  dimensions 
for  the  net,  or  length  dimensions,  we  can  consider  other 
possibilities  for  accurate  designation  of  these  seven  para- 
meters. These  questions,  however,  will  not  be  considered 
here . 
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